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ABSTRACT 

We present a dynamical analysis of the merging galaxy cluster system Abell 2146 using spec¬ 
troscopy obtained with the Gemini Multi-Object Spectrograph on the Gemini North telescope. 
As revealed by the Chandra X-ray Observatory, the system is undergoing a major merger and 
has a gas structure indicative of a recent first core passage. The system presents two large 
shock fronts, making it unique amongst these rare systems. The hot gas structure indicates 
that the merger axis must be close to the plane of the sky and that the two merging clusters are 
relatively close in mass, from the observation of two shock fronts. Using 63 spectroscopically 
determined cluster members, we apply various statistical tests to establish the presence of two 
distinct massive structures. With the caveat that the system has recently undergone a major 
merger, the virial mass estimate is Mvir = x 10 ^“^Mq for the whole system, con¬ 

sistent with the mass determination in a previous study using the Sunyaev-Zel’dovich signal. 
The newly calculated redshift for the system is z = 0.2323. A two-body dynamical model 
gives an angle of 13° — 19° between the merger axis and the plane of the sky, and a timescale 
after first core passage of « 0.24-0.28 Gyr. 

Key words; galaxies: clusters: general; galaxies: clusters: individual: A2146 


1 INTRODUCTION 


Galaxy clusters are formed hierarchically through the mergers of 
smaller clusters and groups. Mergers of galaxy clusters are the most 
energetic events since the big bang, w ith a total kinetic energy that 
can reach 10®^ J (e.g. Sarazin|[2001 i, and major mergers close to 
the plane of the sky are very rare systems but of extreme impor¬ 
tance in cosmology (e.g. |Clowe et al.|2006[|Markevitch & Vikhlinin| 
|2007| l. Most of the mass in galaxy clusters is dark matter, and the 
bulk of baryonic mass is in the form of ionized plasma, comprising 
about 15% of the total mass. During a merger, as the plasma clouds 
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of each cluster pass through each other they are affected by ram 
pressure, causing them to slow down. The galaxies in a cluster are 
effectively collisionless and are affected primarily by tidal interac¬ 
tions during the merger. Dark matter does not have a large cross- 
section for self-interaction (e.g. |Markevitch et al.| 2004{ |Randall et | 
|al.|2008^ , so shortly after collision the plasma clouds are expected 
to lag behind the dark matter and the major concentrations of clus¬ 
ter galaxies (e.g. |Clowe et al.||2006| ). A merger thus results in an 
offset between the dominant luminous component, imaged using 
X-ray telescopes, and the total mass, mapped using gravitational 
lensing. 

Abell 2146 is in the throws of a major merger ( [Russell et al.| 
|2010|l and X-ray maps show a structure similar to the “Bullet Clus- 
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Figure 1. Left: Exposure-corrected Chandra X-ray image for the 0.3-7 keV energy band (North is up and East is to the left). The image has been smoothed 
with a 2D Gaussian of width 1.5 arcsec. Right: HST F435W, F606W, F814W, colour composite image of Abell 2146. 


ter” ( [Markevitch et al.|[2002[ |Clowe et al.||2006| l. Unique among 
merger systems, two clear shock fronts are detected in the Chandra 
X-ray temperature and density maps indicating Mach numbers of 
~ 2.3 for the how shock and « 1.6 for the upstream shock, and cor¬ 
responding velocities in the plane of the sky of 2700jltnn kms“^ 


and 2400 ± 300 km s ^ respectively (Russell et al. 


2012 1 . From the 


shock velocity and the projected distance from the collision site, 
|Russell et al.|p01^ estimate that the merger is observed 0.1 — 0.2 
Gyr after core passage. The detection of the shock fronts is fully 
consistent with the merger being observed soon after first core pas¬ 
sage (e.g. H arazin pool) , but the timescale is approximate due to 
uncertainty in the velocity of the subcluster - likely to be signifi¬ 
cantly lower than the shock velocity (e.g.|Spring el & Farrar|2007l 
|Milosavljevic et al.|2007| l - and since the location of the collision 
site is not precisely known. 


In Fig. 1 we show the Chandra X-ray observatory image from 
|Russell et al.| ( |2012^ and Hubble Space Telescope (HST) optical 
colour composite of Abell 2146 (Program 12871, PI King). Sub¬ 
clusters A and B are labeled on the left, with subcluster A being 
located at the head of the ‘bullet’ in the south-east. The Brightest 
Cluster Galaxy (hereafter BCG) of the SE subcluster (A2146-A), 
is in an unexpected location lagging behind the X-ray cool core 
( [Canning et al.|20f^ . As noted above, in the early stages of a clus¬ 
ter merger the galaxies should be nearly collisionless and follow 
the dark matter distribution, which should lead the collisional X- 
ray gas (e.g. [Clowe et al.|2006[|Markevitch & Vikhlinin|2001^ . 


In addition to the strange projected location of the BCG in 
Abell 2146-A, the system is unusual in that despite it appearing to 
undergo a major galaxy cluster merger it does not have a detected 
radio halo ( [Russell et al.[[20fT) . Studies of samples of radio ha¬ 
los in clusters have shown that they are found almost exclusively 
in morphologically disturbed clusters, indicating that their genera¬ 
tion is associated with acceleration or re-acceleration of energetic 
particles in the highly shocked cluster gas, and that there is a re¬ 
lationship between the radio power of the halo, Pradio, and the X- 
ray luminosity of the cluster, £x ^Liang et al.[[2000[ Giovannini 


[& Feretti][2000| [Buote[[200T] [Brunetti et al.[[2009| [Rossetti et al. 

[201 l[[Cassano et al.|2010b] l. The conditions under which radio ha¬ 
los are generated in clusters, such as the cluster masses, velocities 
and timescales, are still open questions. Merging clusters, such as 
Abell 2146, which do not exhibit these features may be a powerful 
diagnostic to test models of radio halo generation. 


In this paper we present multi-object spectroscopy of the sys¬ 
tem, in order to identify substructure and to provide a comple¬ 
mentary method of probing the dynamical status, timescale of the 
merger, and mass of the system (e.g. [Barrena et al.|2002{ [F errari| 
[et al.|2005][Girardi et al.|1998) . Section 2 will briefly describe the 
optical sample selection, and observations and data reduction will 
be the subject of Section 3. Results from an analysis of the sub¬ 
structure within the Abell 2146 system will be presented in Section 
4 and Section 5 will discuss these results in the context of simple 
dynamical models. Section 6 presents a summary and discussion 
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Figure 2. Target selection for GMOS spectroscopy: Targets were chosen based on colour, magnitude and cluster-centric distance, detected from our Subaru 
Suprime-Cam imaging. Left: Colour-magnitude diagram for galaxies within 3 arcmin of the centre of the system.The dashed line indicates the initial apparent 
magnitude cut which we imposed on our sample. The red sequence was determined by htting a Gaussian in magnitude bins of 18 > Rc, 18 < Rc < 19.5, 
19.5 < Rc < 21 and 21 < Rc. The final 170 galaxies selected from our catalogue by the GEMINI mask making tools are indicated in red and the two BCGs 
are marked by black stars. Right: The 97 objects for which we determined redshifts overlaid on our Subaru Suprime-Cam composite colour image; those 
found to be cluster members are indicated in red and those not in the cluster indicated in white. The circles in each case indicate galaxies which fall on the red 
sequence while the boxes indicate galaxies blueward of the red sequence. 


of our results. A future paper will address the stellar populations of 
the system. 

Throughout this paper we assume a ACDM cosmology with 
present day Hubble parameter Ho = 70kms“^ Mpc“^, matter 
density parameter flm = 0.3, and dark energy density parameter 

= 0.7 (we assume a cosmological constant, equation of state 
parameter w = —1). For this cosmology and at the redshift of 
Abell 2146 {z = 0.2323), an angular size of 1 arcsec corresponds 
to a physical scale of « 3.677kpc. 


2 OPTICAL SAMPLE SELECTION 

Here we briefly present the criteria for selecting the galaxy sample 
for GEMINI GMOS-N multi-object spectrograph (MOS) observa¬ 
tions of Abell 2146. The next section will outline the observations 
and data reduction. 

Targets were detected on our Rc band optical Subaru imag¬ 
ing using SEXTRACTOR ( |Bertin & Amouts||1996| >. We restricted 
our target galaxies to those brighter than 21 apparent magnitude 
and to those within 3 arcmin from the centre of the Abell 2146 
system. We then determined the cluster “red-sequence” by binning 
the galaxies into magnitude bins of 18 > Rc, 18 < Rc < 19.5, 
19.5 < Rc < 21 and 21 < Rc mags. We fit two Gaussians to the 
bimodal population of colour in each magnitude bin; one defining 
the “red sequence” and one defining the “blue cloud”, and took the 
3a boundaries of the Gaussian defining the red sequence to deter¬ 
mine its width in each magnitude bin. We then fit a linear relation 
to the red sequence (black solid lines in Fig. [^. We identify 337 
targets with Rc< 21 mags on the cluster red sequence (see Fig. 2, 
left). 

Many galaxies morphologically associated with the X-ray 
shock fronts and one of the BCGs do not lie on the red sequence. 


For this reason we additionally add 50 targets to our target list (in¬ 
cluding 12 acquisition targets) based on an inspection of the Sub¬ 
aru images. In particular we added the south-east BCG and targets 
close to the BCGs and X-ray shock front which were too blue to be 
included on the red sequence. We also added two possible strong 
lensing features to our target list. Finally we have a sample of 375 
target objects and a further 12 acquisition objects to input to the 
Gemini MOS Mask Preparation Software (GMMPS). 

The two BCGs, galaxies close to the shock front and strong 
lensing features were assigned the highest priority. All other targets 
were prioritized by magnitude then cluster-centric distance and the 
GMMPS was iteratively run in order to fit as many objects as pos¬ 
sible on the masks while requiring at least 3 acquisition objects per 
mask. Finally, we obtained 7 masks totaling spectra on 170 targets 
indicated by red points on the left hand panel of Fig. 2. Of the 170 
targets observed with GEMINI 139 lie on the red-sequence giving 
a red sequence completeness of 67% at Rc< 20.5 mags within 3 
arcmin of the center of the Abell 2146 system. As described in Sec¬ 
tion 3, we were able to obtain redshifts for 97 galaxies, of which 
63 were determined to be cluster members and are tabulated in the 
Appendix (see Appendix and Section 4). The 97 objects are indi¬ 
cated on the right panel of Fig. 2 with the cluster members in red. 
Objects indicated by circles are those determined to be on the red 
sequence while those indicated by boxes are found blueward of the 
red sequence. 


3 OBSERVATIONS AND DATA REDUCTION 

Spectroscopic observations of Abell 2146 were carried out with the 
Gemini Multi-object Spectrograph (GMOS) at the Gemini North 
Telescope on Mauna Kea, Hawaii. They were made over the course 
of four nights between 2012 April 19 and May 1 (Program ID: 
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Figure 3. One dimensional extracted spectra of the BCGs in A2146-A (left) and A2146-B (right). Spectra are shown at the observed wavelengths. Various 
emission and absorption features are labelled. 


GN-2012A-Q-47, PI Canning). We used the GMOS-MOS instru¬ 
ment equipped with the B600-l-_G5307 grating. This configuration 
yielded a dispersion of O.SOA/pixel. The 5.5 arcmin field of view 
of the CMOS detector allowed observations to be made in a single 
pointing. 

Using slit widths of 1 arcsec and slit lengths of 5 arcsec, be¬ 
tween 18-31 objects were able to be targeted in each mask. Due to 
the densely packed nature of the field, 7 masks were used to obtain 
spectra of a total of 170 science targets. Three 1200 second expo¬ 
sures were taken with each mask in order to achieve a S/N~ 3 per 
spectral pixel, at the redshift of the 4000A break for the faintest 
galaxies in our target list (determined using GEMINI CMOS expo¬ 
sure time calculator). Wavelength calibration was performed using 
Helium-Argon lamps. 

The data reduction was performed using the Gemini package 
in IRAiQ see Fig. 3 for examples of the resulting spectra. All sci¬ 
ence exposures, comparison lamps, and flats were overscan/bias 
subtracted and trimmed. The resulting reduced spectra were then 
processed to remove cosmic rays, wavelength calibrated, extracted 
to ID and combined. 

The cross-correlation method was used to calculate radial ve¬ 
locities through the RVSACj^ package in IRAF. A catalogue of 23 
different templates from the RVSAO library ( [Kurtz & Mink|199^ , 
including stars, quasars, and emission and absorption line galaxy 
template spectra was used to determine the radial velocities. In or¬ 
der to determine the templates to use in the final analysis we relied 
on a combination of visual examination and the R-value returned 
by RVSAO. Each extracted ID target spectrum was visually exam¬ 
ined and a redshift identified by at least two members of the team, 
with a third examining spectra where disagreements were found. 
For 54/170 objects no redshift was able to be determined with a 
visual inspection due to a low S/N. For a further 19/170 targets 

^ IRAF is distributed by the National Optical Astronomy Observatories, 
which are operated by the Association of Universities for Research in 
Astronomy, Inc., under cooperative agreement with the National Science 
Foundation. 

^ RVSAO was developed at the Smithsonian Astrophysical Observatory 
Telescope Data Center. 


the visual examinations yielded a redshift which was not consistent 
with the radial velocity measurements for any of the 23 RVSAO 
template spectra. In both these cases the targets were not used in 
the remainder of the analysis in this paper. In 73/170 targets the 
visual examination yielded a redshift within ±6,000km s“^ of the 
RVSAO radial velocity with the highest R-value. In these cases the 
radial velocity corresponding to the highest R-value template was 
used. In the remainder (24/170), the highest R-value template did 
not correspond to the visually determined cz. However, a template 
with R-value smaller than the best value did have a radial velocity 
within ±5,000 km s“^ of the visual velocity. If A i? < 5 then 
we accepted the fit with a lower R after further examination of the 
template spectra. In most of these 24 targets poor sky subtraction 
was the reason for poor template fitting. However, the HST imag¬ 
ing has shown many lensed features are apparent in the Abell 2146 
system and galaxy density is high. Some poor redshifts may be due 
to blending between two of more sources. The spectra of the BCGs 
in A2146-A and A2146-B are shown in Fig. 3. 

In order to minimize errors associated with cross-correlation 
and to correct for artifacts in the spectra (such as cosmic rays), three 
separate observations of each object were made and then the spectra 
were co-added together (e.g. jGallazzi et al.|2008[ . This signal aver¬ 
aging helps to increase the signal to noise of the spectrum and can 
help to remove certain undesirable features that are only contained 
in one of the individual spectra. The radial velocities and errors 
were then compared to both the co-added and individual spectra. 

4 A2146 STRUCTURE 

4.1 Cluster Member Determination 

The target radial velocities were used to determine cluster mem¬ 
bership: first, we made a ±10,000 kms“^ cut at the previously 
calculated cluster redshift of z = 0.2343 ( jStruble & Rood|1999| (, 
as was done by jFerrari et al.| ( [^05| (. Cluster members were then de¬ 
termined by a 2.5(7 clip of the remaining objects. This leaves a total 
of 63 galaxies out of the 97 objects for which we could determine a 
redshift. Fig. 2 shows the 63 objects identified as cluster members 
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Figure 4. Top: The radial velocity histogram for all objects where a reli¬ 
able redshift was determined and the targets are not stars. The cluster radial 
velocities peak around 70,000 km s“^, however, there is no evidence in 
our data that the broad tails in the distribution correspond to further clus¬ 
tering of galaxies. Bottom: Radial velocity distribution of cluster members. 
Cluster members in A2146-A are shown in blue, and A2146-B in red. The 
solid lines indicate a Gaussian profile at the biweight mean. The radial ve¬ 
locities of the two BCGs are denoted by stars. Subcluster membership was 
determined using the KMM algorithm ^Ashman et al.|1994}. 


circled in red. The foreground and background objects are identi¬ 
fied in white. The velocity distribution for the 97 targets (which are 
not stars) is shown in the top panel of Fig. 4 while the bottom panel 
shows only the 63 cluster members after a clipping. 


4.2 Subclustering 

The velocity distribution of the Abell 2146 system is shown in 
Fig. El We apply a Lilliefors test to determine whether the velocity 
distribution of the 63 cluster member galaxies is statistically dif¬ 
ferent from a Gaussian distribution. We calculate a p-value of 0.23 
indicating that the null hypothesis, that the velocity distribution of 
the system as a whole is Gaussian distributed, cannot be rejected. 
Flowever, for a plane of sky merger this is not a sufficient condition 
to exclude the presence of subclustering. 


66.382 
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Q 
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Figure 5. Projected density map of cluster members in 2D. The relative 
density contour map is shown in position space. The subcluster A2146-A 
is on the bottom left (blue circles) and A2146-B (red circles) is on the top 
right with the BCGs in both subclusters represented with a black point. In 
the figure North is up and East is to the left. 
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To investigate the possibility of subclustering further we apply 
the Dressler-Schectman test ( [Dressier & Shectman|1988l > (hence¬ 
forth DS test) and also carry out an analysis of the 2D projected 
density of cluster members. 

First, we apply the DS test which provides a useful statistic 
to determine the likelihood of subclustering. The DS test compares 
the kinematics of a sample of size Nnn to the kinematics of the 
whole system. For each galaxy in the field, the distances are com¬ 
puted between every pair in the system and then the Nnn nearest 
neighbours are chosen for analysis. The statistic is given by 

^2 ^ “ vf + {cTlocal “ Cr„)^] , (1) 

LZ y 

where (t„ is the velocity dispersion of the whole group and aiocai is 
the velocity dispersion of the sample. These values are summed 
to obtain 


N 

A = y] <5*, (2) 

i 

with ^ > 1 being an indicator of substructure present in the field. 
Using Nnn = \/NTotai we found a. ^ — 1.08. The marginal result 
makes it hard to determine if there is substructure in the system (see 
also [Pinkney et al.]jl996^ for the results of applying the DS test to 
simulated galaxy clusters). 

To test if the calculated A was significantly different from 
random, we randomly shuffled the velocities of the objects while 
holding the positions fixed. Repeating this process 10,000 times, 
we found that the fraction of random velocities that gave a larger A 
was P = 0.18, where a value of P less than 0.1 would be expected 
for a non-random distribution. Thus, we conclude that the DS test 
results are not strong enough to indicate substructure in the system. 
We return to this result in the Discussion. 

Second, we analyse the 2D projected density of all of the 
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cluster members. Cluster members identified by their redshift were 
plotted in RA and DEC space with overlaying contours determined 
using a Kernel Density Estimator (KDE) and a Gaussian kernel. 
The KDE utilized a Silverman bandwidth. Since Gaussian kernels 
are less sensitive to bimodal distributions, this test will be biased 
towards strong substructure detections. The results plotted in Fig. 5 
clearly show that two density peaks are present in the data. This 
confirms the existence of at least two subclusters and is consistent 
with the Chandra X-ray data which shows a bi-model merger mor¬ 
phology with both forward and reverse shocks detected in the X- 
ray analysis and along the same axis as our two subcluster clumps 
( [Russell et al.|2012) . The KDE analysis allows us to conclude that 
there is substructure present in the system. 

Having confirmed the presence of substructure we wish to 
identify the cluster members of particular subclusters in order to 
characterise properties such as the mass of the substructures, the 
timescale since core passage of the merger and the presence of 
any offsets between the galaxy and hot gas populations. In order 
to assign the cluster members to particular subclusters we apply the 
KMM algorithm ( [Ashman et al.|199^ . This determines the mixture 
of multi-dimensional Gaussians which are best fit to a set of data 
and the probability of each data point belonging to each Gaussian. 

The clear X-ray shock fronts and the small line-of-sight ve¬ 
locity offsets of the two BCGs indicates strongly that our merger 
is likely to be in the plane of the sky. As such, the distribution of 
velocities from each subcluster is expected to have a significant 
overlap (as seen in Fig|^. Including such similar velocity informa¬ 
tion in the KMM algorithm can lead to non-physical results in dis¬ 
tinguishing between the two subcluster components. The analysis 
was carried out with position data only and yielded results which 
correspond well with the projected galaxy density peaks seen in 
Fig. 0 The KMM algorithm classified 35 galaxies as belonging to 
A2146-A and 28 belonging to A2146-B. These subcluster mem¬ 
bers are marked in Fig. (with blue corresponding to A2146-A 
and red corresponding to A2146-B), and subcluster membership is 
also noted in Table 1 in the appendix. The velocity distribution of 
the galaxies identified in the two subclusters is shown in Fig.l^ The 
KMM algorithm is sensitive to outliers ( [Ashman et al.[1994t which 
is of particular importance for a cluster sample such as ours where 
the data is incomplete. We therefore bootstrap the sample (e.g. Wa¬ 
ters et al. 2008) and re-run the KMM test 5,000 times to determine 
which galaxies to assign to which subcluster and the centroid of our 
two subclusters. The centroids and their 68 per cent uncertainties 
are 239.062tg °gt 66.349tg;°“ (15:56:14.9lg;®, 66:20:56.4l^J) 
and 239.005tg;|;g|, 66313t.^ °o°o 2 (15:56:1.21};^ 66:22:22.8i}J) 
for subclusters A and B respectively. 

Applying a Lilliefors test again to the individual subclusters 
yields a p-value of 0.83 for A2146-A and a p-value of 0.059 for 
A2146-B. This indicates that the velocity distribution for A2146-A 
can be adequately described by a Gaussian distribution. A Gaus¬ 
sian distribution for the velocities of subcluster A2146-B would be 
rejected at the 94 (< 2cr) per cent level. There is evidence from 
a weak gravitational leasing mass reconstruction of HST data that 
A2146-B is rather more extended that A2146-A, possibly due to an 
intrinsically more elongated mass distribution (King et al. 2015). 
The numbers of galaxies for which we have spectroscopic data is 
insufficient to decide between these possibilities. 


4.3 Velocity Dispersions of the subclusters 

The biweight estimator jBeers et al.|1990^ was applied to the clus¬ 
ter members for location and scale giving a mean apparent cluster 


redshift of 2 = 0.2323 (radial velocity of 69, 640 km s“^). The ve¬ 
locity dispersions were calculated from the median absolute devia¬ 
tion (MAD) (e.g. [Ferrari et al.|2005^ with errors estimated through 
a bootstrap technique. The velocity dispersions of the two subclus¬ 
ters are 1130l}2okm s“^ for the A2146-A subcluster and TOOl}®® 
kms“^ for the A2146-B subcluster, with a 68% confidence inter¬ 
val (Cl). It is important to note that the distribution of radial veloci¬ 
ties is under-sampled and given that A2146 has recently undergone 
a merger, the distributions are possibly non-Gaussian. So the cal¬ 
culated velocity dispersions characterise the distribution, with the 
caveat that the sample may be non-Gaussian. 

Application of the biweight estimator and bootstrap error ap¬ 
proximations were carried out using the astropy library in the 
Python programming languag^ In a bootstrap technique, a 68% 
Cl is found for the MAD. This approach draws many samples with 
replacement from the data available, estimates the MAD from each 
sample, then rank orders the means to estimate the 16 and 84 per¬ 
centile values for 68% Cl. Unlike assuming normal distributions 
to calculate 68% Cl, the results calculated from the bootstrap are 
robust and can compensate for a data distribution that is far from 
normal. 

The velocity distribution of cluster members is shown in the 
histogram in Fig. 4. The cluster members in A2146-A are shown 
in blue, and the cluster members in A2146-B are shown in red. 
Gaussian profiles for each subcluster are overlaid and centered at 
the mean determined from a biweight estimator ( [Beers et al.[1990[ l. 


5 A2146 DYNAMICS 
5.1 Virial Mass 


A common way to estimate the mass of a galaxy cluster is through 
the virial theorem, which assumes dynamical stability (e.g.[Araya-[ 
[Melo et al.|2009) , and therefore relates the time averaged kinetic 
energy to the potential energy as 


1 ^ 

= (3) 

i 

We calculate the virial masses below and tabulate them in Tablejl] 
along with the SZ and X-ray derived masses from the literature (see 
Section 6). However, we caution that the Abell 2146 system is not 
dynamically relaxed and discuss the potential biases of the mass 
measurements in Section The virial theorem relates the gravita¬ 
tional potential energy to the kinetic energy (relative velocities) of 
the cluster members and the mass is calculated via the following 
relation 

= (4) 

where a„ is the velocity dispersion calculated in Section [43] Rh 
is the mean harmonic radius given by 


1 

Rh 


2{N-2)\ 

iV! 


E 



(5) 


where rij is the separation vector between the and galax¬ 
ies and N is the total number of objects. The benefit of using a 
mean harmonic radius is that it more accurately measures the ef¬ 
fective radius of the gravitational potential of the cluster members 


^ (http://docs.astropy.org/en/stable /1 
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( |Araya-Melo et aL]|2009[ l, important when considering the virial 
state of a system. This radius also has the benefit of being indepen¬ 
dent of the cluster centre allowing it to reflect internal structure of 
the system, putting extra weight on objects that are closer together 
( |Araya-Melo et al.|2009^ . The virial radius for the entire system is 
related to the mean harmonic radius by r^ir ~ 2 Rh ~ 0.64Mpc 
(e.g. |Girardi et al.|1998[[Merchan & Zandivarez|2005^ . 

Taking the values calculated for the velocity dispersions of 
each of the subclusters along with the harmonic radius for each, 
68% Cl virial mass estimates are found to be Mvir = 6.OI2 g 
10^^ Mq for A2146-A and ? x 10^"^ Mq for A2146- 

B (we note here that subcluster B was found not to be well de¬ 
scribed by a Gaussian distribution), giving a mass ratio of 2.4. As 
an estimate for the total mass, we take the sum of the two estimates, 
since the assumption that the individual clusters are virialized is 
likely more appropriate than assuming that the system as a whole is 
virialized. This yields Mvir = x Note that the er¬ 

ror bars do not reflect systematic errors which can change the mass 
estimate, even when the system is not in the throws of a merger; 
this simple dynamical estimate assumes a virialized spherical mass 
distribution with isotropic velocity dispersion. 

The large errors in the mass estimates are due to the large er¬ 
rors inherent to the velocity dispersions. Since the virial mass is 
proportional to the square of the velocity dispersion, the error on 
the mass estimates is quite large. The 68% Cl level corresponds to 
a lc7 uncertainty. For a discussion of the impact of systematic and 
model dependent errors on the mass estimates, see Section]^ 

Since the estimates of mass from the kinematics of galaxies 
make several simplifying assumptions, we consider now an alterna¬ 
tive estimate of the mass of each subcluster derived from a scaling 
relation between velocity dispersion and mass. |Evrard et al.]j2008^ 
used N-body simulations to obtain a prediction for the relationship 
between mass and velocity dispersion of dark matter halos. The ad¬ 
vantage of this type of approach is that it includes more complex 
dynamical behaviour. Although their simulations were carried out 
using only dark matter, since galaxies and dark matter particles can 
both be treated as collisionless particles moving in a cluster po¬ 
tential, they are expected to have approximately the same velocity 
dispersion ( [Okamoto & Hab e|1999[ see a lso discussion in |Evrard] 
|et al.|200^ . Using the Evrard et al.|p008| l relation 

M 2 ooE{z) = 9.358 X x , (6) 

with the parameter E{z) = 1.264 at the cluster redshift and for the 
cosmology we work in, and determining the mass estimate for each 
cluster gives a total mass of M 200 = ® x 1Q2^Mq. 

The observed non-Gaussianity of A2146-B in conjunction 
with the system not being dynamically relaxed will contribute to an 
overestimation of the velocity dispersion for each subcluster com¬ 
ponent. The |Evrard et al.| ( |2008| l scaling relation does not account 
for these and therefore it is likely that the mass will be over esti¬ 
mated as well. 

Specifically for merging clusters, there have been studies us¬ 
ing simulations to investigate the accuracy of various ways of esti¬ 
mating mass jPinkney et al.|199^ Zuhone et al.|2009) or focusing 
on specific merger syste ms (e.g. [Zuhone_elj^ |2009^ . For clusters 
with a mass ratio of 3:1, [Pinkney et al. 1 1996l show that the mass 
inferred from a simple application of the virial theorem, and when 
sub-clustering is ignored in the analysis, can be biased by as high 
as a factor of two. The maximum bias arises when the clusters are 
merging with an axis along the line of sight, and is somewhat less 
in mergers with an axis close to the plane of the sky, though still 
up to a factor of 1.5. The hydrodynamic simulations of [Zuhon"^ 



Figure 6. Projected angle from the plane of the sky, a, plotted against the 
difference in radial velocities of the the subclusters, V%. The black vertical 
line gives the solutions for Vu = 620kms“^ and the red for Vr = 850 it 
60kms“^. The red shaded area around 3% represents the errors in radial 
velocity from the cross-coirelation. The dotted curve represents the limits of 
the bound solutions as calculated from Eq. 8. The gray region represents the 
area where only unbound solutions would exist. The intercept of the vertical 
lines gives solutions corresponding to a = 13° — 19°, V = 2600 — 
2800kms-l,andil = 0.70 - 0.72Mpc. 


|et al. 


, _ 

( 1996^ . Zuhone et al.| ( |2009| l also support accounting for separate 

subclusters when determining mass estimates from the velocities of 
galaxies in a merger system. 


consistent with the earlier work of|Pinkney et al.| 


5.2 Two-Body Dynamics 

We undertake a dynamical analysis of the system using a two-body 
model, making the assumption that our system consists of two dom¬ 
inant subclusters, which is reasonable based on the projected den¬ 
sity map (see Fig.l^. We use the methods adopted by Beers et al 
( [Beers et al.|1982| i7to which the reader is referred for further de¬ 
tails, and also follow the formalism presented in e.g. Girardi et ay 
l [1998^ , [Barrena et ^j2002[ l and [Ferrari et al.[ ( [2005| . Taking the 
total mass of the system derived from the virial mass estimates and 
the difference in radial velocities of the two subclusters, we calcu¬ 
lated the angle between the merger axis and the plane of the sky, 
and the age of the merger. Note that this method assumes point 
masses on radial orbits with no net rotation. The subcluster halos 
will also most likely be overlapping in a recent merger. Making 
these assumptions in this analysis will add systematic and model 
dependent uncertainties to the results obtained. Therefore the solu¬ 
tions will be approximations to the dynamics of A2146. The cou¬ 
pled equations, which describe the time evolution of the system, are 
closed by setting t = to = 3.4 x 10^^ s, the age of the universe at 
the Abell 2146 system redshift 2 ; = 0.2323. 

The parametric solutions to the equations of motion are given 
as follows: 

-R = ^(1 - cos(x)); (7) 
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V Rm / 1 - cos(x) 

where R is the spatial separation, V is the relative velocity between 
the two BCGs, Rm is the separation of the two subclusters at max¬ 
imum expansion, M is the total virial mass of the system, and x is 
the so-called developmental angle. 

The limits of the bound solutions can be found by considering 
the gravitational binding energy in the following relation 

Vr/Zp ^ 2G'MTotai sin^ (a) cos(a) (10) 

where Vr is the difference in radial velocities of the two subclusters 
and Rp is the projected separation of the cores of the two subclus¬ 
ters. 

Since the velocity distributions of the subclusters, A2146-A 
and A2146-B, do not precisely follow a Gaussian distribution, the 
difference in the mean radial velocity of the subclusters from a 
Gaussian fit (from Fig.|^ will not yield an accurate value of Vr. 
We therefore opt to use the difference in radial velocities of the two 
BCGs in each of the subclusters as a determination of Vr. 

The spectra of the two BCGs are shown in Section 3 (Fig[^. 
The BCG in A2146-B has a strong absorption spectrum with well 
defined features, as is typical of BCGs. The spectrum of the BCG 
in A2146-A, however, shows many emission lines. The three most 
prominent features are Hp and two [OIII] lines. Zooming in on 
this spectrum (Fig|^ shows the characteristic absorption spectrum 
expected for an elliptical galaxy with additional strong emission 
features. 

While most of the absorption features are blended with emis¬ 
sion features, making it difficult to accurately determine redshift, 
the Na D (rest wavelength of 5895 A) feature is relatively un¬ 
blended and yields a velocity of ~69,612 km s“^. This is a lower 
redshift than determined through cross correlation methods for the 
BCG which rely on the emission lines (69,897 km s“^). 

|et al.|pOT2l l show that much of the emission line gas in the A2146- 
A BCG is offset in a plume from the centre of the galaxy which 
may be the cause of the discrepancy in the radial velocities of the 
emission and absorption lines. 

Calculating Vr as the difference between the absorption line 
spectrum of A2146-B and the emission spectrum of A2146-A gives 
a value of Vr = 850 ± 60 kms“^. Using instead the Na D ab¬ 
sorption feature yields a value of Vr = 620 kms“^. Nominally, it 
would be advantageous to take the velocities of the absorption lines 
of the galaxies, however, as the BCG in subcluster A is also offset 
from its predicted location - leading, rather than trailing, the cool 
core, which may indicate that the BCG is not in this case a good 
tracer of the center of the mass distribution - we opt to consider 
the uncertainty as the relatively broad range in velocities from the 
emission and absorption spectra. 

We therefore take Vr to be in the range 620 km s“^ to 850 ± 
60kms“^ in the two-body analysis. Rp is taken as the projected 
separation of the two BCGS. At the cluster redshift this corresponds 
to 0.69 Mpc. 

The radial velocity and projected radius are related to the sys¬ 
tem parameters as follows: 


Canning 


Vr = Vsin{a)\Rp = Rcos{a). (11) 

Solving equations IZ10 and 1^ iteratively yields the solutions 
shown in Fig. 6. Plotting Vr shows the solutions for our system. 
While there are solutions of a close to 80°, the presence of the ob¬ 
served shocks in the X-ray analysis ([Russell et al.|2010^ and the es¬ 
timated value of 17° ( [Canning et al.|2012 1 , suggests that the merger 
is close to the plane of the sky. This allows us to take the range of 


A2146-A BCG 



Figure 7. Zoomed in view of the spectrum for the BCG in A2146-A. Both 
absorption and emission features are clearly visible and the radial velocity 
of the emission features are offset by Ri -1-230 km s“^ from the BCG. Both 
absorption and emission features are labeled. 


solutions for a = 13° to 19°, which are very close to the values 
derived by [Canning et al.| ( |2012[ . The results from the two-body 
model analysis are V = 2600 — 2800km s“^ for the relative ve¬ 
locity, and R = 0.70 — 0.72 Mpc for the separation. Assuming 
a collision site of « R/2, half way between the two subclusters, 
yields a merger age of 0.24 — 0.28 Gyr. 


6 DISCUSSION 

X-ray observations of Abell 2146 have revealed the system has two 
clear shock fronts and a gas structure similar to the Bullet cluster 
( [Markevitch et al.[2002) showing it must be in the throes of a major 
merger. The detection of two clear shock fronts also indicate that 
the merger should be in an early stage. However, the Abell 2146 
system presents a couple of challenges to that hypothesis. 

The first is that the A2146-A subcluster BCG is located behind 
the bow shock and in the wake of the ram pressure stripped cool 
core rather than leading it (see[R ussell et al.|2010[ and [Canning et[ 
[al.|2012[ l. As mentioned in Sectionshortly after a galaxy cluster 
collision the intra-cluster medium is expected to lag behind the bulk 
of the mass distribution traced by the dark matter and the major 
concentrations of cluster galaxies (e.g. [Clowe et al.[2006^ . 

One possibility for the unusual location of the BCG could be 
that the cluster merger is not in its early stages (e.g. [Hallman*^ 
[Markevitch[2004[ (. In the later stages of a galaxy cluster merger the 
collisionless matter distributions slow due to gravity allowing the 
gas to catch up and be accelerated in a ‘gravitational sling shot’ 
beyond the dominant mass distribution. Other possibilities are that 
the system may not be adequately represented by a two-body col¬ 
lision with the cluster containing complex substructure, or that the 
merger is somewhat off-axis. 

Examining the redshifts of our target galaxies we find 63 
galaxies are likely members of the A2146 system. An analysis of 
the projected density peaks of the galaxies (see Section [4!^ shows 
that the Abell 2146 system is dominated by two subclusters. Our 
data do not provide evidence for other more complicated structures 
such as other merging groups or subclusters. We use the KMM al- 
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gorithm to assign subcluster membership to the two subclusters and 
find 35 subcluster member galaxies in A2146-A and 28 in A2146- 
B. Application of the DS test to determine substructure was in¬ 
conclusive. However, we emphasise that in [ Dressier & Shectman| 
( [1988^ and in jPinkney et al.| ( ri996| l, the sensitivity of the DS test to 
substructure in merger systems is shown to be diminished when the 
subclusters are of similar mass and when their merger axis is close 
to the plane of the sky, as is likely the case in Abell 2146. 

Given that the system appears to be fairly well described by 
two dominant subclusters we undertake a simple two-body analy¬ 
sis of the subcluster velocities and positions (section [5!2| and Fig.|^ 
to determine the properties of the merger. As the velocity distribu¬ 
tions of the subclusters are not both well described by Gaussian 
distributions we choose to use the BCGs as a tracer of the centre 
of the subcluster mass distributions. We note here that in the case 
of subcluster A2146-A this may be incorrect since, as mentioned 
above, the BCG is offset from the expected location of the centre 
of the mass distribution; whilst the centroid of the distribution of 
red sequence galaxies is approximately coincident with the X-ray 
cool core, the BCG appears to lag behind the cool core. If a colli¬ 
sion has offset the BCG we may therefore expect the BCG velocity 
not to accurately trace the subcluster velocity centroid. The A2146- 
A BCG not only shows absorption features but also strong emis¬ 
sion features in the spectrum which have a velocity offset of ~230 
km s“^. This is significantly larger than the errors in the velocities 
from cross-correlation (~30 km s“^). The centroid of the subclus¬ 
ter velocity distribution falls in between the emission and absorp¬ 
tion values (A2146-A; BCG emission 69,897 km s“^, BCG absorp¬ 
tion 69,612 km s"\ Centroid (KMM) 69,878 km s”!; A2146-B; 
BCG absorption 68,827 km s“^. Centroid (KMM) 69,710 km s“^; 
WYFFOS spectra published in jCanning et al. [201 2 [ determined sim¬ 
ilar velocities of 69,753 and 68,940 km s“^, from single Gaussian 
fits to the emission and absorption spectra, for BCG A and B re¬ 
spectively). We therefore choose to use the full range of velocities 
from the absorption and emission lines in the A2146-A BCG in the 
two body analysis to allow a conservative estimate of the merger 
axis angle and timescale. 

Despite these uncertainties in cluster velocity centroid, the 
merger age, for bound solutions close to the plane of the sky, is 
still predicted to be young (0.24—0.28 Gyr), similar to that deter¬ 
mined for the Bullet cluster jBarrena et al.|20()^ . The subclusters 
are therefore likely being observed shortly after first core passage. 
The unusual BCG position is therefore unlikely to be due to a grav¬ 
itational ‘slingshot’ of the gas past the galaxies (see e.g. [Hallman [ 
[& Markevitch|2004l (. 

The merger axis is determined to be very close to the plane of 
the sky with a — 13° to 19°. The additional a solution (top of Fig. 
^ that would correspond to a line-of-sight merger is unlikely due to 
Russell et al.H2010|[2012[ > observing both forward and reverse X- 
ray shock fronts. Observations of merger shocks are rare in X-ray 
observations as they are easily diluted by projection effects, they 
are therefore most prevalent in systems which are close to the plane 
of the sky. We therefore conclude that the unusual BCG location is 
not due to projection effects in a line-of-sight merger. 

These findings for the merger timescale and merger axis are 
also consistent with the presence of the two large shock fronts seen 
in Chandra X-ray observations jRussell et al.|2012l l, since later in 
the merger the shock fronts would not be visible, and for a larger 
angle with respect to the plane of the sky their sharp edges would 
be smeared by projection effects. 

Although the two-body model is useful for estimating the dy¬ 
namical history of merger systems, the assumptions that the clus- 
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^500 [^q] 

Figure 8. Distribution of clusters in the Pi .4 — M 500 diagram (Fig. 3 
from Cassano et al. (2013)). The blue circles were changed to diamonds for 
clarity and additional lines are added corresponding to dynamical and X- 
ray analysis. The various different symbols indicate halos belonging to the 
extended Giant Metrewave Radio Telescope (GMRT) Radio Halo Suivey 
jCassano et al.|20T3) , or EGRHS, (blue diamonds); halos from the litera¬ 
ture (black circles); halos with very steep spectra (USSRH, green asterisks); 
A1995 and Bullet cluster (blue stars); cool core clusters belonging to the 
EGRHS (magenta arrows). Best-fit relations to giant RHs only (black lines) 
and to all RHs (including USSRH, green dashed lines) are reported. The 
circled blue arrow is A2146 with an M 500 and Pi .4 reported by Cassano 
et al. (2013). The dashed line shows A2146 and it’s derived M500 from 
dynamical analysis. This corresponds to a radio halo power upper limit of 
Pi .4 Si 4 X lO^^W Hz“^ for a mass of M500 = 4.3 X The 

horizontal purple line marks the calculated 3(7 upper limit for Abell 2146 
from |Russell et aQj2011) . 

ters are point masses with constant mass over time, start to break 
down when they begin to merge. The analysis also does not in¬ 
clude a change in relative velocity of the clusters with time. In the 
context of the Bullet Cluster, [Nusser[ (2008[ > discusses several lim¬ 
itations of the two-bodel model, for example how dynamical fric¬ 
tion and tidal stripping are not included in this type of analysis. 
Both of these effects act to reduce the relative speeds of the clus¬ 
ter components. Although hydrodynamic simulations are the most 
accurate way to model the dynamics of merger systems, there are 
prohibitive in terms of the complexity of setting up the simulations, 
and the computational expense. [Dawson| ( [2013^ presents a novel 
approach where the computational requirements fall somewhere 
between the two-body model and hydrodynamic simulations. By 
comparison with full hydrodynamic simulations of the Bullet Clus¬ 
ter, this method is accurate to about 10% in determining parameters 
describing the clusters and their merger trajectory. Required inputs 
to this dynamical modeling approach include spectroscopic red- 
shifts for cluster members, arising from this work, and additional 
mass estimates available from the weak lensing analysis (King et 
al. 2015). Thus application of the [Dawson| ( [2013^ method to Abell 
2146 will be the topic of future work. 

The second challenge for Abell 2146 is its apparent lack of 
radio emission, neither radio relics, which would be expected to be 
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associated with the shock fronts, nor a radio halo have been seen 
( [Russell et al.||20fT| l. Large-scale diffuse radio emission in clus¬ 
ters has only been found in merging systems and as such it could 
be a useful diagnostic of the statistics of dynamically relaxed and 
merging systems. However, the mechanisms which form this dif¬ 
fuse radio emission and the conditions required for their formation 
are not yet known. 

Diffuse radio emission on Mpc-scales in clusters requires local 
particle acceleration where either a pre-existing electron population 
is re-accelerated by merger turbulence (e.g. [Brunetti et al.||2001| 
|Petrosian|2001| ( or emitting electrons are continuously injected by 
collisions between cosmic ray protons and thermal protons (e.g. 
jPennisonl 19801 for reviews of extended radio emission in clusters 
see jPeretti & Giovannini||2008l jFerrari et al.||200^ jPeretti et al.j 

|2012| l. The recent discovery of extremely steep radio halo spectra 
(e.g. [Brunetti et al.|2009[ jPallacasa et al.|2Q0 ^ together with the 
non-detection of nearby clusters at 7 -ray wavelengths (e.g.jBrunettij 
jet al.|201^|Ackerman et al.|2010t favours turbulent re-acceleration 
models over the ‘secondary’ continuous injection models. 

Correlations between the synchrotron power in radio halos, 
fradio, and the galaxy cluster X-ray luminosity. Lx and temper¬ 
ature, Tx, have been found, which also show a bi-modailty for 
merging and non-merging galaxy clusters (e.gJLiang^££^ 2000 


Kemnper & Sarazin|2001[[Bacchi et al.|2003l|Cassano et al. 2006 


Brunetti et al. |2007> . Clusters with detected radio halos have Pradio 
correlated with Lx and show a disturbed morphology consistent 
with merging where there is no radio halo detected, and with only 
upper limits on radio luminosity, the clusters generally appear more 
relaxed and lie roughly an order of magnitude below this corre¬ 
lation jBrunetti et al.|200^ [Rossetti et al.|20irt . However, some 
clusters, such as Abell 2146, have been found which confound this 
bi-modality in the X-ray (e.g. Russell et al.[[20lT] [Giovannini &| 


[Peretti[2000[[Cassano et al.|2010b[ l. 


Russell et al. 


([2011[| showed the 3cr radio halo upper limit is a 


factor ~5 times lower than predicted given its X-ray luminosity. X- 
ray luminosity is boosted in a merger and this factor is greater for 
young mergers jRicker & Sarazin|2001^ so the confirmed young 
age of the Abell 2146 system (being only 0.24 — 0.28 Gyr old) 
and the near equal masses of the subclusters indicates its X-ray 
luminosity is likely to be enhanced by a factor of ~2-5. A boost this 
significant could place Abell 2146 back on the relation for merging 
clusters. However, the age since core passage is similar to other 
systems in which radio halos have been detected. 

[Brunetti et al.[ ( (T012[ ( examine a sample of radio halo clus¬ 
ters from the GMRT Radio Halo Survey and find that cluster syn¬ 
chrotron emission is suppressed on timescales of only a few hun¬ 
dred Myrs. The simulations of [Donnert et al.[ ( [2013[ l, albeit for a 
1:7 mass ratio merger, suggest that the radio halo luminosity should 
be increasing rapidly during infall and through core passage, with 
the X-ray and radio luminosity, within a few hundred Myr, being 
boosted by factors of 2 and 30, respectively. They find the radio 
emission starts to decay roughly 0.5 Gyr after core passage. While 
timescale may be relevant in boosting the X-ray luminosity and 
therefore in where Abell 2146 should sit on the Pradio-Lx relation, 
these and other simulations suggest that the timescale since core 
passage for Abell 2146 is not too short or too long for the genera¬ 
tion of a radio halo. 

Cluster mass is thought to be a key component of radio halo 
generation in galaxy clusters ^Cassano & Brunetti[[2005] l. [Basu[ 
\2Q\2) examine the radio-SZ correlation for clusters on the basis 
that the integrated SZ signal is a more robust indicator of the cluster 
mass regardless of dynamical state than the X-ray luminosity (e.g. 


[Mod et al.|2005[[Nagai et al.|2007^ . A correlation is found for the 
radio power and integrated SZ signal, Ysz for the sample of clus¬ 
ters and no bi-modality is detected between those clusters which 
are relaxed and those which are disturbed. The authors suggest the 
bi-modality with X-ray properties may be due to a bias in the radio 
quiet. X-ray selected clusters, being selected towards lower mass 
cool core clusters. However, they note that their sample size for ra¬ 
dio quiet systems is small and so a weaker bi-modality might still 
exist. jCassano et al.[2013) l find the converse, that a bi-modailty ex¬ 
ists, in both the Pradio-Lx and Pradio-Psz relations with all relaxed 
systems having radio halo upper limits which sit >2a below the 
relation. However, they suggest mass is also important and that the 
bi-modality only exists for ‘SZ-luminous’ Ysz > 6 x 10“® Mpc^, 
or alternatively massive, M 500 > 5.5 x Mq galaxy clusters. 

The scaling relation of [Evrard et al.[ ( [^08[ l, while robust, does 
not account for the observed subclustering or dynamical state of 
Abell 2146. While our individual subcluster masses similarly suf¬ 
fer from the dynamics of the merger the total mass as the sum of 
the subcluster masses rather than assuming one virialised system 
should yield value more in line with previous estimates. We there¬ 
fore use a total mass of Mvir = x 10^^ Mq giving a mass 

ratio of approximately 2.4:1 (A2146-A to A2146-B). The large un¬ 
certainties in the mass estimates preclude determining which sub¬ 
cluster is in fact the more massive component with a high degree of 
confidence. However, the system appears to have undergone a re¬ 
cent relatively equal-mass merger. The X-ray analysis found that 
the bow shock and upstream shock are of comparable strength, 
which is consistent with a close to equal mass merger, however 
a 1:4 mass ratio with A2146-B being more massive than A2146-A 
was favored [Russell et al. 1 20K)|(. For high er mass ratio systems, 
such as the Bullet cluster ( [Clowe et al.[2006t with a ratio of ~ 10:1, 
the perturbation to the gravitational potential of the merger is much 
smaller and therefore will produce a much weaker upstream shock. 
Our dynamical mass ratio estimate is consistent with a weak tens¬ 
ing analysis of the system (King et al. 2015). 

The NFW (Navarro, Frenk & White 1996) r 2 oo value of the 
system is approximately equal to the virial radius, or to twice 
the harmonic radius calculated in Section 5.1, yielding a value of 
r 200 = 0.64 Mpc. Assuming an NFW mass profile and a halo con¬ 
centration of 4 typical of massive galaxy clusters (see collation 
from King & Mead 2011), an estimate of rsoo « 0.42 Mpc was 
calculated, corresponding to M 500 = Q x 10^"^ M©. 

To take into account the overestimate in mass due to the sys¬ 
tem having undergone a recent merger, a surface pressure term 
(SPT) can be subtracted from the viral mass to get a more accu¬ 
rate estimate of the true mass of the system ( [The & White|19^ . 
Assuming the velocities are isotropic, the SPT correction can be 
taken as SPT = 0.2 x Mvir ( [Girardi et al.|1998) . The reduction 
of the mass due to the SPT correction would result in a lower mass 
estimate of M 200 = 6.8 x 10^"^ Mq and M 500 = 4.3 x 10^'* M©. 

For Abell 2146, SZ observations ( [Rodri g uez-Gonzalvez et[ 
[ar][201 l| l yielded a total system mass inside r 2 oo of M 2 oo,sz = 
4.1 ± 0.5 X M© (or 5.61 q'y x Mq using the same 

cosmology as this paper i.e. h = Ho/70 kms“^Mpc“^). This 
is a factor of ~1.5 smaller than our dynamical virial mass esti¬ 
mate for the total system mass not including a SPT correction, or 
roughly 1.2 times smaller with the SPT correction, bringing the 
calculated M 200 within the errors of the SZ mass estimate. The 
M 500 and M 200 values calculated from the SZ (taking h = Ho/70 
kms“^Mpc~^), our dynamics analysis, and the X-ray analysis of 
[Russell et al.[p01^ are compiled in Table 1. 

When correcting M 500 for the SPT, and using the direct scal- 
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Table 1. Summary of Mass Estimates for Abell 2146. All masses are in 
units of Mq . 


Mass 

M500 

M 20 Q 

Dynamics 

5.4I3 3 X IQi"' 


X 1014 

sz 

4.3^°;® X 1014 


X 1014 

X-ray 

7tl X 1014 




ingof M500 — Pi A, shown in Fig. 7 (Fig. 3 from |Cassano et ^ 
2013 1 , we find an expected radio halo power of Pi ,4 « 4 x 10^^ 


Watt FIz~ , marked with the dashed black horizontal line in Fig. 7. 
The solid purple horizontal line marks the calculated 3cr upper limit 
for A2146 from the GMRT radio observations by [Russell et al.| 
( |201 1^ . This observational upper limit is therefore very similar to 
the radio power predicted by the observed M500 — Pi. 4 correlation 
for merging systems, using a dynamical mass estimate. 

A lower M 500 value could be the leading factor as to why 
no extended radio emission is detected, and this is accommodated 
within the errors on the dynamical mass, particularly if we include 
the SPT. As the radio observing sensitivity increases in the future, 
faint extended radio emission may well be detected and the possi¬ 
bility of a mass threshold for radio halos in merging clusters will be 
clarified. Furthermore, the prospects for assessing the dependence 
of radio halo production on other factors will improve with studies 
of large samples of major cluster mergers. 


7 CONCLUSIONS 


We present a dynamical analysis of the galaxy cluster merger sys¬ 
tem Abell 2146 using multi-object spectroscopic observations from 
the Gemini North Telescope. We employ subshucture tests such as 
the Dressler-Schectman test ( [Dressier & Shectman [198^ , analy¬ 
sis of an adaptive kernel density map ( Ptsam^^W , and applica¬ 
tion of the KMM algorithm j Ashman et al.|[1994^ using the po¬ 
sition information of the galaxies to conclude that the system is 
dominated by two main subshuctures, and is likely well described 
by a two-body collision. One cluster is located in the south-east 
(A2146-A) and the other in the north-west (A2146-B). An average 
redshift of z = 0.2323 was calculated with velocity dispersions of 
llSOj*! 320 km s“^ and 76 OI 340 kms“^ respectively for the A2146- 
A and A2146-B subclusters. The main results are: 


(i) From the 63 spectroscopically confirmed A2146 system 

members, we estimate a total virial mass of Mvir = x 10 ^^ 

Mq , without any correction applied for the ~20 per cent boost in 
mass due to the system having undergone a recent merger (e.g.|The| 
[& White|1986l >. 

(ii) The two subclusters have a close to 2.4:1 mass ratio indicat¬ 
ing that this is a major merger. 

(iii) A two-body analysis of the system yielded an angle be¬ 
tween the merger axis and plane of sky of approximately 13° — 19° 
and a time since core passage of the merger of approximately 
0.24 - 0.28 Gyr. 

These findings for the dynamics are consistent with the obser¬ 
vation of large shock fronts on Chandra X-ray images (Russell et 
al. 2010; 2012). We conclude from our analysis that the unusual lo¬ 
cation of the BCG, behind, not in front of the shock front and X-ray 
cool core, is not due to the merger being in a late stage of its evolu- 
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tion with the gas experiencing a gravitational sling-shot beyond the 
main mass distribution. 

Despite the lack of an observed radio halo, simulations of our 
calculated timescale and near equal mass ratio of the merger would 
seem to favor the generation of a radio halo (e.g. [Donnert et al.[ 
[2013[[Brunetti et al.|2012[ (. However, the final mass of the system, 
and therefore necessarily the progenitor masses of the two subclus¬ 
ters, is fairly low, compared with other similar cluster mergers. [Cas-[ 
[sano et aT] ( [2013[ l suggest that a bi-modality in radio halo proper¬ 
ties for merging and non-merging systems may require a total mass 
threshold of about M500 > 5.5 x Mq. Uncertainties on the 
mass estimates are very large, however, Abell 2146 would lie close 
to this threshold. Deeper radio data on this and other merging sys¬ 
tems that appear to lack radio halos would be critical for testing 
whether a threshold in mass exists. Magnetohydrodynamic simu¬ 
lations of the system, constrained by the wealth of observational 
data, will further our understanding of the distribution of luminous 
and dark matter in the system, as well as the radio emission. For the 
Bullet Cluster, recent work by [Lage & Farrar[p014[ ( illustrates that 
this approach provides much insight into merger systems. Appli¬ 
cation of the [Dawson| \20 1 3) method for dynamical modeling will 
provide a more accurate description of the geometry of the merger 
than provided by the two body model, also greatly narrowing down 
the parameter space for MHD simulations. 

Conflicting scenarios have been proposed where merger 
events may be a catalyst to spur star formation - for example shock 
fronts moving through the ICM and triggering further starbursts 
( [Mihos & HernquTst|[1996[ [Bekki & Couch[[2003^ - or to quench 
star formation via ram-pressure stripping (Fujita et al. 1999). De¬ 
tailed analysis of the [O II], [O III], Ha, and TT/j spectral lines is 
currently underway (Lee et al. in prep), in order to estimate the star 
formation rates in the cluster galaxies as discussed in [Ferrari et al.j 
( [2005^ . 
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Table 1: Data for 63 spectroscopically confirmed cluster member galaxies in the 
field of Abell 2146. The radial velocities are calculated from RVSAO cross¬ 
correlation methods applied to our Gemini CMOS spectroscopic observations 
(PI Canning). 


ID 

RA 

(degrees) 

Dec 

(degrees) 

CZ 

(km s“^) 

Acz 

(km s~^) 

Subcluster 

(A/B) 

Type 

(Emission/Absorption) 

5 

238.9647 

66.37298 

69484 

24 

B 

Ab 

13 

238.9836 

66.36637 

68014 

130 

B 

Ab 

14 

238.9839 

66.37659 

69128 

47 

B 

Ab 

16 

238.9883 

66.38978 

67482 

35 

B 

Ab 

21 

238.9973 

66.37249 

70186 

16 

B 

Ab 

22 

238.9987 

66.36988 

67705 

10 

B 

Em 

23 

239.0007 

66.37329 

68849 

24 

B (BCG) 

Ab 

24 

239.0027 

66.36994 

71123 

37 

B 

Ab 

35 

239.0064 

66.36687 

68517 

17 

B 

Ab 

36 

239.0085 

66.3711 

71163 

31 

B 

Ab 

39 

239.013 

66.3703 

71325 

21 

B 

Ab 

43 

239.0155 

66.37942 

68398 

19 

B 

Ab 

47 

239.0175 

66.377 

70649 

20 

B 

Ab 

48 

239.0196 

66.35606 

70406 

31 

B 

Ab 

50 

239.022 

66.3901 

69195 

58 

B 

Ab 

54 

239.0246 

66.36791 

71930 

49 

B 

Ab 

55 

239.025 

66.3782 

68802 

21 

B 

Em 

57 

239.0291 

66.36815 

71835 

37 

B 

Ab 

60 

239.031 

66.3541 

70504 

78 

A 

Ab 

61 

239.034 

66.349 

69249 

78 

A 

Ab 

63 

239.034 

66.3764 

69432 

32 

B 

Ab 

67 

239.0385 

66.3488 

71634 

14 

A 

Ab 

68 

239.0385 

66.3294 

68242 

33 

A 

Ab 

71 

239.0403 

66.36081 

72876 

16 

A 

Ab 

76 

239.0445 

66.3462 

71054 

81 

A 

Em 

79 

239.0495 

66.36106 

69924 

27 

A 

Ab 

80 

239.0505 

66.3496 

71365 

25 

A 

Ab 

81 

239.0505 

66.3534 

69371 

43 

A 

Ab 

83 

239.052 

66.3567 

72327 

41 

A 

Em 

87 

239.0535 

66.3701 

69925 

38 

A 

Ab 

90 

239.056 

66.33381 

68368 

54 

A 

Ab 

92 

239.0565 

66.3565 

67810 

37 

A 

Ab 

94 

239.058 

66.3482 

69897 

26 

A (BCG) 

Em 

96 

239.0601 

66.34456 

70033 

22 

A 

Ab 

105 

239.067 

66.3461 

69497 

24 

A 

Ab 

108 

239.0691 

66.34259 

68075 

21 

A 

Ab 

109 

239.0694 

66.34006 

71013 

52 

A 

Ab 

110 

239.07 

66.3467 

68929 

16 

A 

Ab 

113 

239.0762 

66.3713 

68790 

28 

A 

Ab 

115 

239.0768 

66.36112 

69552 

31 

A 

Ab 

116 

239.0795 

66.3415 

71610 

18 

A 

Ab 

117 

239.0799 

66.34054 

69059 

32 

A 

Ab 

120 

239.0827 

66.34473 

69121 

39 

A 

Ab 

122 

239.0835 

66.3427 

72489 

42 

A 

Em 

123 

239.0835 

66.3564 

71683 

35 

A 

Ab 

128 

239.0955 

66.3482 

70673 

31 

A 

Ab 

135 

239.1241 

66.35526 

68352 

18 

A 

Ab 

138 

238.989 

66.3766 

68624 

30 

B 

Ab 

139 

239.0091 

66.36326 

68960 

22 

B 

Ab 

142 

238.9803 

66.36338 

69323 

23 

B 

Ab 

143 

238.9758 

66.36388 

66669 

25 

B 

Ab 

147 

239.0616 

66.33845 

69860 

14 

A 

Ab 

149 

239.0802 

66.33684 

67965 

25 

A 

Ab 

150 

239.085 

66.343 

70321 

16 

A 

Ab 
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Table 1 - Continued 


ID 

RA 

(degrees) 

Dec 

(degrees) 

CZ 

(km s“^) 

/S.CZ 

(km s~^) 

Subcluster 

(A/B) 

Type 

(Emission/Absorption) 

152 

239.0406 

66.34549 

67638 

18 

A 

Ab 

153 

239.058 

66.356 

70031 

22 

A 

Ab 

154 

239.04 

66.3502 

70282 

17 

A 

Ab 

156 

239.0025 

66.3748 

69179 

27 

B 

Ab 

157 

239.016 

66.3708 

68956 

27 

B 

Ab 

164 

239.0625 

66.3677 

71160 

39 

A 

Em 

165 

239.0184 

66.36444 

69028 

31 

B 

Ab 

168 

239.0051 

66.39083 

70413 

35 

B 

Em 

170 

238.959 

66.3771 

69239 

23 

B 

Ab 
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